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Abstract. We present a photometric method for identifying stars, galaxies and 
quasars in multi-color surveys and estimating multi-color redshifts for the extra- 
galactic objects. We use a library of > 65000 color templates for comparison with 
observed objects. The method was originally developed for the Calar Alto Deep 
Imaging Survey (CADIS), but is now used in a variety of survey projects. We 
checked its performance by spectroscopy of CADIS objects, where it provides high 
reliability (6 mistakes among 151 objects with R < 24), especially for the quasar 
selection, and redshifts accurate within ct^ ~ 0.03 for galaxies and ~ 0.1 for 
quasars. For an optimization of future surveys, a few model surveys are compared, 
which use the same amount of telescope time but different filter sets. In summary, 
medium-band surveys perform superior to broad-band surveys although they collect 
less photons. A full account of this work is already in print ([11,12]). 

1 Introduction 

The value of a multi-color survey is usually not simply limited by the depth 
of its imaging but by the depth of its classification used to select the objects 
of interest from the survey data. Here, we investigate what factors the depth 
of a successful classification depends on. Obviously, it depends in some way 
on the filter set used as well as on the depth of the imaging, but the optimum 
choices might vary according to the goal of the survey. 

If a survey aims at identifying only one type of object with characteristic 
colors, a tailored filter set can be designed. E. g., when looking for U-band 
dropouts [9] , the ugr filter set is a very good choice. The performance of such 
a dropout survey depends mostly on the depth reached in the U-band, so the 
photon flux detection limit in U is the key figure. Also, number count studies 
are limited by the completeness limit in the filter of concern. 

If all kinds of different objects are the subject of the survey, a filter set 
covering the entire spectral range accessible to the survey instrument is a 
natural choice. The SLOAN Digital Sky Survey [13] is presently the most 
ambitious project to provide a broad-band ugriz color database, on which the 
astronomical community might perform a large number of virtual surveys. 

Of course, it is possible to split the bands of a filter set even further into 
a larger number of medium-band filters. Surveys like CADIS with typically 
10 to 20 filters are sampling the visual spectrum with a resolution almost 
comparable to that of low resolution imaging spectroscopy. Intuitively, one 



2 Christian Wolf et al. 



might expect that such surveys are very expensive in terms of telescope time 
while their particular value might not be obvious. This paper will show, that 
in fact medium-band surveys are a powerful and economically competitive 
alternative to broad-band surveys. 

2 Classification and Redshift Estimation in CADIS 

CADIS fostered the development of a scheme for spectral classification, that 

distinguishes stars, galaxies, quasars and strange objects. Simultaneously, it 
assigns multi-color redshifts to extragalactic objects, an idea dating back to 
1962 [1]. Essentially, a likelihood-based algorithm compares the observed col- 
ors with a library of template colors. The probability of each single template 
to generate the colors of the observed object is calculated and a probability 
density function is obtained for each object class, dependent on the redshift 
and the spectral energy distributions (SED) contained in the library. From 
these functions, relative likelihoods for class memberships and estimates for 
redshift, SED types and their errors can all be derived (see [11] for all details). 

The color template libraries arc assembled from observed spectra by syn- 
thetic photometry performed on our filterset. As an input we used the stellar 
library by Pickles [8], the galaxy template spectra by Kinney et al. [5] and 
the QSO template by Francis et al. [3]. From this, we generated regular grids 
of quasar templates ranging in redshift within < z < 6 and having various 
continuum slopes and emission line equivalent widths. The quasar spectra are 
attenuated bluewards of the Lyman-a line by a rcdshift-dcpendent through- 
put function modelling the Lyman forest. Also, a grid of galaxy templates 
has been generated for < < 2, and contains various spectral types from 
old populations to starbursts. 

Using this method, the object database of CADIS was classified and the 
result was checked by spectroscopy for 151 objects with R < 24. An object 
with i? = 24 is typically detected at a 5-a level in most of the up to 17 filters. 
At i? < 22 the classification was almost free of mistakes with two mistakes 
among 103 objects (two Seyfert-l-QSOs were classified as normal galaxies). 
In the fainter regime of i? > 22, a morphological criterion was applied to all 
extended objects classifing them directly as galaxies and four classification 
mistakes were finally found among 48 objects. 

Multi-color redshifts were accurate to ~ 0.03 for 90% of the galaxies at 
i? < 24 while the remaining 10% are estimated significantly wrong. Among 
the quasars, 50% were accurate to <Jz w 0.1 with another 50% of complete 
mistakes. At the present time, the observed quasar color spectra are affected 
by intrinsic variability during the four years of CADIS observations needed to 
collect the data for all the bands. Once this effect is corrected for we expect 
improvements for the multi-color redshifts of the quasars (see [12] for details 
on the CADIS results). 



Classification and Redshift Estimation 3 



Table 1. Filters and lO-a-magnitude limits for the three survey setups (A, B, C) 
compared with Montc-Carlo simulations. The I-band filter is a long wavelength 
passband filter with a cut-on wavelength roughly at 780 nm. Its far-red sensitiv- 
ity limit is given by the dropping quantum efficiency of the CCDs. All filters are 
installed in the Wildficld Imager at the 2.2m-MPG/ESO telescope at La Silla ob- 
servatory. Setup A has actually been used for the deep WFI survey COMBO-17. 



Acen/fwhm (nm) 
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23.5 


23.5 


24.1 


456/99 
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25.0 


25.0 


25.6 
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24.5 


24.5 


25.1 


652/162 
850/150* 
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24.5 


24.5 


25.1 
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23.0 


23.0 


23.6 
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23.58 
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23.1 
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23.0 






696/20 




22.8 


23.18 




753/18 




22.7 






815/20 




22.6 


22.98 




856/14 




22.5 






914/27 




22.4 


22.78 





3 Model Surveys in Comparison 



Initially, it should be natural to assume that surveys with different filter sets 
show quite a different performance in terms of classification and redshift es- 
timation. If a survey aims for objects with very particular spectra, the filter 
set can certainly be tailored to this purpose. If the objects of interest span 
a whole range of spectral characteristics, it is not trivial to guess via ana- 
lytic thinking which filter set performs best. Here, we present a comparison 
of three fundamentally different filter sets and show their performance for 
classification and redshift estimation (sec [12] for all details). 

The three model surveys spend the same total amount of exposure time on 
different filter sets, using the same instrument, telescope and observing site. 
We chose the Wild Field Imager (WFI) at the 2.2-m-MPG/ESO-tclcscopc 
on La Silla as a testing ground, because it provides a unique, extensive set 
of filters ranging from several broad bands to a few dozen medium bands to 
choose from. Furthermore, the WFI is a designated survey instrument which 
is extensively used by the astronomical community. 
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The three modeUed surveys, here called setup "A" , "B" and "C" , each 
spend ISOksec of exposure time distributed on the following filters (see also 
Tab.|l|): 

Setup A spends 50 ksec on the five broad-band filters of the WFI (UBVRI) 
and 100 ksec on twelve medium-band filters. Using ESO's exposure time cal- 
culator V2.3.1 for the WFI, we related exposure times to limiting magnitudes 
assuming a seeing of 1.4 arcsec, an airmass of 1.2, point source photometry 
and a night sky illuminated by a moon three days old. The exposure times 
are distributed such, that a quasar with a power-law continuum f^, = 1/°' and 
a spectral index of a = —0.6 is observed with a uniform signal-to-noise ratio 
in all medium bands. As a result, the twelve medium bands each deliver a 
10-fT detection of an i? = 23.0-quasar. Setup B spends 50 ksec on the same 
broad bands but concentrates the 100 ksec for medium-band work on only six 
filters reaching a uniform lO-cr detection of a i? = 23.38-quasar then. Setup 
C finally spends all 150 ksec on the broad-band filters and omits the medium 
bands entirely. 

The resulting performance was evaluated in terms of the fraction of suc- 
cessfully classified objects as a function of class, magnitude and survey setup 
and also in terms of the achieved accuracy of the multi-color redshifts. As 
shown in Fig.|l|, all surveys perform rather equal in terms of class recovery. 
At i? = 21 ... 23, the surveys work almost perfect and then degrade quickly 
with magnitude until becoming useless at i? > 24. Obviously, the setups with 
medium-band filters included (A, B) collect much less photons than the pure 
broad-band setup (C), but deliver a comparable classification performance. 
Therefore, the smaller number of detected photons is compensated by the 
increased information content per photon. 

In contrast, the accuracy of the multi-color redshifts shows very significant 
differences between the setups (Fig.||). In the bright regime at R < 22, the 
redshift resolution is proportional to the wavelength resolution of the filter 
set, setting the medium- band surveys far ahead the broad-band survey. The 
information potential of medium-band filters, which do not average the object 
spectra as broadly as broad-band filters do, can fully be exploited and cause 
a clear advantage for the total performance of medium-band surveys. Only in 
the faint regime at i? > 24, where the redshift estimation is almost useless, 
the differences between the setups vanish. 

These results on the performance are average figures for large samples. In 
this context it is important that pure broad-band surveys contain a few class 
ambiguities in multi-color space. After dropping the u-band from the SDSS 
filter set (which provides only a weak detection of faint objects anyway), an 
M2V star and a galaxy of type Sb at z = 0.5 can not be distinguished by 
color, even if the relative calibration of the filters is accurate to 3%. Similarly, 
a QSO at z = 3 with a — —0.8 displays the same colors as a G2V star. In 
contrast, the two medium-band filter sets can easily distinguish these pairs 
of objects to rather faint levels. 
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Fig. 1. Fraction Q of simulated objects wiiidi are correctly classified in the three 
different setups (solid line = setup A, grey line = setup B, dashed line = setup 
C) over the R-band magnitude in a f 50 ksec-survey at the WFI@2.2 (see Tab. f for 
details). Basically, all setups are comparably successful (taken from [11]). 
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Fig. 2. Variance of redshift estimation error {Az = Zmc — z) among simulated 
objects in the three different setups (lines as in Fig. 1). Setups A and B provide 
the highest redshift resolution. Early type galaxies work better due to their higher 
continuum contrast at the 4000 A-break. 



4 Information Content of Multi-Color Surveys 

The rather equal performance of the different survey setups can be under- 
stood by analytically considering an infomation content for the surveys, given 
by the product of the following factors: 

• number of filters n 

• exposure per filter t cx n^^ 

• filter width AX oc n^^ 

• information per detected photon $ cx AX~^ oc n 

The total information content yields Itot oc ri° =const. This result is only 
valid under the ideal assumption that photon noise was the only source of 
error and the color libraries mimic true nature perfectly. But in a real survey 
there will be calibration errors limiting the achievable best accuracy to 3% 
(equivalent to a 30-a detection) and also cosmic variance will cause observed 
spectra to deviate in a scattered manner from the library templates. In the 
bright regime, where performance is limited by these systematic effects and 
not by photon statistics, many noisier datapoints provide more information 
than fewer precise but inaccurate ones. So, under realistic circumstances the 
medium-band survey setups are likely to be more powerful than the pure 
broad-band setup even when being limited by equal telescope time. 
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5 Ongoing Applications 

We currently work on a number of multi-color surveys, where filters and ex- 
posure times were chosen to match some primary survey strategy. Although, 
none of these might have been optimal choices in terms of a general classifi- 
cation, we used or intend to use our approach to extract class and redshift 
data on the objects contained. These surveys are: 

1. The Calar Alto Deep Imaging Survey (CADIS): Four broad-band (BRJK') 
and twelve medium-band filters have mostly been chosen to match the 
needs of the emission line survey in CADIS, while some of them fill in gaps 
in the spectral coverage. The multi-color part of CADIS has been used to 
study the evolution of the galaxy luminosity function at z = 0.1 . . . 1.3, 
to search for quasars at all visible rcdshifts and to use the observed faint 
stellar population to check models of the Galactic structure and the stel- 
lar luminosity function [4,6,7,10,12]. 

2. A lensing study of the galaxy cluster Abell 1689: Two broad-band and 
seven medium-band filters have been chosen to separate well between the 
cluster galaxies at z w 0.19 and the background population. The galaxy 
luminosity function in the background of the cluster is compared to a 
control field taken from CADIS, and the cluster mass is estimated from 
weak lensing effects on the apparent luminosities [2]. 

3. A WFI survey for Classifying Objects by Medium-Band Observations with 
17 filters (COMBO-17): The filters (setup A from Tab. 1) are chosen 
to provide a selection function and a redshift accuracy for quasars and 
galaxies, which is as independent of redshift as feasible. The data will 
be used to study the faint end of the quasar luminosity function at all 
accessible redshifts z > 1, galaxy-quasar correlation at z < 1.3, as well 
as weak lensing effects in the cluster group Abell 901/2 and in the open 
galaxy field. 

4. The Sloan Digital Sky Survey (SDSS): Five broad filters are used to span 
the entire range of available CCD sensitivity. We intend to apply our 
classification to get multi-color redshifts, to search for quasars and to 
separate stars from compact galaxies, where morphology data are not 
sufficient. 

6 Summary 

We presented an innovative method that performs a multi-color classification 
and redshift estimation of astronomical objects in a unified approach. The 
method is essentially based on templates and evaluates the statistical consis- 
tency of a given measurement with a database of spectral knowledge, serving 
as a second, very crucial input to the algorithm. Spectroscopic confirmation 
within CADIS showed that the method works well and is of high practical 
relevance. 
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A simulation of a few model surveys with fundamentally different filter 
sets revealed that in practice surveys with many medium-band filters are 
superior to pure broad-band surveys for the purpose of classification and 
redshift estimation. Even when being limited by telescope time, the medium- 
band filter sets achieve at least the same working depth as the broad-band 
survey while providing much improved redshift resolution. 

We conclude, that this method should be very suitable for many survey 
applications, which usually require only low spectral resolution and finite 
accuracy in the derivation of physical parameters, but aim for large samples 
to feed statistical studies and to search for rare and unusual objects. Of 
course, if a 100% sure confirmation on the nature of an individual object is 
needed, or if you aim for high resolution studies, the method gives only a 
preselection of candidates. 
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